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ABSTRACT 

We report results from a spectral and timing analysis of M82 X-l, one of the brightest 
known ultraluminous X-ray sources. Data from a new 105 ks XMM-Newton observation 
of M82 X-l, performed in April 2004, and of archival RossiXTE observations are 
presented. A very soft thermal component is present in the XMM spectrum. Although 
it is not possible to rule out a residual contamination from the host galaxy, modelling it 
with a standard accretion disk would imply a black hole mass of « lO 3 Af . An emission 
\ line was also detected at an energy typical for fluorescent Fe emission. The power 

— . i density spectrum of the XMM observation shows a variable QPO at frequency of 11 3 

*-£h ■ mHz with properties similar to that discovered bv IStrohmaver fe Mushotzkvl lf2003'). 

Th e QPO was also found in 7 archiv al RXT E observations, that i nclud e those analyzed 
Q ■ bv IStrohmaver fe Mushotzkvl l)2003|) and iFiorito fe Titarchukl l)2004[) . A comparison 

V- 1 | of the properties of this QPO with those of the various types of QPOs observed in 

c/3 , Galactic black hole candidates strongly suggests an association with the type-C, low 

frequency QPOs. Scaling the frequency inversely to the black hole mass, the observed 
QPO frequency range (from 50 to 166 mHz) would yield a black hole mass anywhere 
in the interval few tens to 1000 Mq. 
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1 INTRODUCTION 

First rev ealed by the Einstein Observatory (see e.g. 
lFabbianolll989l) . point-like, off-nuclear X-ray sources with 
luminosities significantly in excess of the Eddington limit 
for one solar mass h ave been discovered in a large number of 
nearby galaxies (e.g.lColbert fe PtaJj2002)TColbert fe Milled 
l2004l:ISwartz et all2004lLiu fe BregmadbuOSft . The nature 
of most of these ultraluminous X-ray sources (ULXs) re- 
mains unclear. In several cases a variability on timescales of 
months/years has been detected, a possible signature of an 
accreting compact object. Were these sources X-ray binaries 
in the host galaxy, and assuming Eddington-limited accre- 
tion, masses in the range 100-1000 Mq are inferred from 
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the observed flux. The central object would be then an in- 
termediate mass black hole (IMBH). Mass estimates based 
on the Eddington limit are, however, questionable: these 
sources need not to be spherically symmetric, nor station- 
ary. It has been proposed that many of the ULX properties 
can b e explained assum ing that they do not emit isotropi- 
cally dKing et alJbOOfl or are dominate d by emission from 
a relativistic jet fe.g. llCaaret et aT]l200&t) . In this case, they 
may harbor stellar mass BHs and may be similar to Galactic 
black-hole binaries. 

The available XMM-Newton and Chandra observa- 
tions show that, for sources with sufficiently good statis- 
tics, the best fit of the X-ray spectrum is often obtained 
with a two-component model: an absorbed multicolor disk 
(MCD)+power-law. Some ULXs have MCD temperatures 
kTMCD ~ 200 — 300 eV, much lower than in Galactic 
BH binaries. Since Tmcd or MgJ/ 4 , the low tempera- 
ture of the soft component has been taken as further evi- 
dence for the presence of an IMBH of ~ 100-1000 Mq in 
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tiller et al.ll2003l 


2004 ICroDDer et alJl2004 


IZamDieri et all 


2004 iKone et al 


2005). A few sources have 



a temperature above 1 keV, higher than that of Galactic BH 
candidates, but still in the range of Galactic microquasars. 

Optical counterparts of UL Xs have been s tudie d by 
a number of a uthors (see e.g. IZampieri et al] 120041 and 
iLiu et all 120051 and references therein). Several ULXs are 
associated with optical emission nebulae of a few hundred 
parsecs in diameter, whose prope rties are consistent wit h 
those of X-ray irradiated nebulae ijPakull fc Mirionl f2002 ) . 

Another approach to study the nature of ULXs is 
through time variability. The analysis of the aperiodic vari- 
ability in the X-ray flux of X-ray binaries is a powerful 
tool to study the properties of the inner regions of the 
accretion disk aro und compact objects (for a review see 
Ivan der Klisl EoOsT ) . In particular, Quasi-Periodic Oscilla- 
tions (QPOs) provide well-defined frequencies, which can 
be linked to specific time scales in the disk. QPOs can be 
broadly divided into three classes: (a) QPOs at very low fre- 
quencies (<0.02 Hz), probably associat ed to oscillations and 
instabilities in the accre tion disk (see iMorean et al.l Il997t 
iBelloni et aT]ll997t l200Cft : (b) Low-Frequency (LF) QPOs, 
with typical frequencies between 0.1 and 10 Hz, probably 
con nected to similar oscillations in neutron star systems (see 
e.g.lBelloni et alJl2002l: IRemillard et al.ll2002aL Ivan der Khl 
120051: ICasella. Belloni fc Stellal 120051) . over whose origin 
there is no consensus; in Black Hole Candidates (BHCs) 
3 main different types of LF QPOs have been identified 
llCasella. Belloni fc Stellall2005l and references therein); (c) 
"hecto-Hertz" QPOs, with a typical frequenc y of 100-300 
Hz, in two cases observed to appear in pairs l)Strohmaveil 
l2001allbh . It is currently unclear whethe r these QPOs show 
a constant frequency fo r each source fsee lHoman et, all200lt 
IRemillard et all2002bl) . and whether they d o appear in pairs 
obeyin g particular frequency ratios (see IRemillar d et alJ 
2002b). However, since they identify the highest frequen- 
cies observed in these systems, they are the best candidates 
for association with, e.g., the keplerian frequency at the in- 
nermost stable orbit, or relativistic precession frequencies. 

Whatever their physical nature, as they originate in the 
inner regions of accretion disks around black holes, these 
features are expected to be produced also in ULXs. How- 
ever, if ULXs contain IMBHs of 100-1000 M , the fre- 
quencies involved are much smaller. The first and, to date, 
only ULX where a QPO has b een discovered is M82 X-l 
iStrohmaver fc Mushotzkvl2003t) . The QPO has a frequency 
of 54.4 mHz and a FWHM of 11.4 mHz, leading to a quality 
value Q — v/Av ~ 5. The total fracti onal rms of the QPO 
in the 2-10 keV band is 8.4%. Recently. iFiorito fc Titarchukl 
(2004) reported the identification of another QPO at 106 
mHz in the power spectrum of M82 X-l from RossiXTE 
data, arguing that it may be a harmonic of the QPO at 54 
mHz. 

Here we report the results from a joint spectral and tim- 
ing analysis of a 105 ks XMM-Newton observation of M82 
X-l performed in April 2004 and of archival RossiXTE ob- 
servations of the same field. The new XMM data are part of 
a pointing to M82, originally requested to study the enrich- 
ment of the galactic interstellar medium (Ranalli et al., in 
preparation) . The plan of the paper is the following. In S0wc 
present the spectral analysis of t he new XMM observation 
and a re-analysis of the data bv IStrohmaver fc Mushotzkvl 



(2003). reports the long term light curve of M82 X-l. 
In the timing analysis of the new XMM observation and 
other archival RossiXTE observations is presented. Finally, 
in §5\ we discuss our results. 



2 SPECTRAL ANALYSIS 

2.1 XMM-Newton observation of April 2004 

XMM-Newton observed the galaxy M82 on April 21, 2004 
(Observation ID 0206080101, PI: P. Ranalli) for a total ex- 
posure time of 105 ks. The three XMM EPIC cameras (pn, 
MOS1 and MOS2) operated in Full Frame mode with the 
medium filter. Photon lists, data screening, region selection 
and events extraction were performed with standard soft- 
ware (XMM-SAS v 5.4.1 and 6.0.0). Events lists were di- 
rectly extracted from the observation data files using the 
tasks epproc and emproc for the EPIC pn and MOS data, 
respectively. The last part of the exposure was affected by 
high background radiation. Standard high flares filtering (to- 
tal off-source count rate above 10 keV less than 1 count s^ 1 
for pn and 0.35 count s _1 for MOS) leaves 76/72 ks of good 
time intervals for the pn/MOS cameras. Slightly different 
good time intervals were used for the spectral and timing 
analysis, as specified below. 

Source counts of M82 X-l were extracted from a circu- 
lar region of radius 13", centered on the coordinates RA= 
09:55:50.2, DEC=+69:40:47. A previous observation with 
Chandra showed several sources i n or nearby the position of 
M82 X-l jMatsumoto et al .1200 ll) . To avoid as much as pos- 
sible contaminations we subtracted the contribution of some 
of them by excluding circular regions of 5" radius centered 
on the Chandra positions (see Figure^ . The nearest sources 
could not be eliminated in this way. Another strong source 
of background conta mination, especially at low e nergie s, is 
the host galaxy itself. IStrohmaver fc Mushot zkv ( 2005J) did 
not attempt to subtract it and limited their analysis to en- 
ergies ^ 2 keV. Here we try to perform spectral fits of the 
XMM data in the [0.8-10] keV range and therefore perform 
subtraction of the host galaxy diffuse emission. The back- 
ground was extracted from an annulus of inner and outer 
radii of 18 and 35", respecti vely, eliminating the cont ribu- 
tion of sources 1, 8 and 9 of Ma tsumoto et al] i200 1]) in a 
way similar to that described above. The background sub- 
tracted count rate is ~1.0 count s _1 for pn and ~0.35 count 
s" 1 for each MOS. 

Spectra were grouped to require at least 150/80 counts 
per bin for pn/MOS and analyzed with XSPEC v. 11. 2.0. 
Standard interstellar absorption (wabs) was taken into ac- 
count and an overall normalization constant was used to 
minimize residual detector calibration uncertainties. Results 
of the spectral fits with various models are reported in Ta- 
bic Q The spectrum is rather composite, so that fits with 
single component models do not give satisfactory results. In 
particular, the residual contribution from the host galaxy 
at soft energies (< 2 keV) needs to be a d equat ely mod- 
elled. In a different context, IStevens et al] i2003T) tried to 
fit the nuclear emission of M82 with two MEKAL (optically 
thin, thermally emi tting plasma) compo nents plus a power- 
law (PL). Following IStevens et all (12003^1 . we tried to fit the 
residual galactic contamination with a dual MEKAL model 
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Figure 1. XMM EPIC MOS1 exposure of M82, taken on April 21, 2004. Left: full image of the galaxy in the [0.2-10.0] keV band. Right: 
l'xl' region centered on M82 X-l in t he energy range [2.0-10. 0] keV. Contour plots show the position of all the sources resolved in the 
Chandra observation of January 2000 iMatsumoto et al 1l200ll) . The big circle (13") shows the region used for source count extraction, 
while the two smaller (5") circular regions are used to subtract the contribution of sources 2 and 3 of Matsumoto et al. (2001). 
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Figure 2. XMM EPIC pn spectrum of M82 X-l in April 2004 
(Obs. ID 0206080101) with the best fitting model {heavy solid 
line) and components: PL {gray solid line) + MEKALs {thick 
dotted lines) + MCD {thick solid line). 



with different absorbing columns. Using this model for the 
diffuse emission plus an absorbed PL for the ULX gives a 
reasonable fit (x 2 = 1.31 for 846 d.o.f.). This can be fur- 
ther improved adding a third MEKAL component, needed 
mainly to fit a residual excess in emission around 6.6 keV 
(x 2 = 1-24 for 843 d.o.f.). This line, attributed to fluorescent 
Kq Fe emission, can also be fitted adding a gaussian to the 
PL component. 

Adding an additional soft component to the PL results 
in a further improvement of the fit. The best fit of the joint 
pn and MOS spectra of M82 X-l is obtained adding to the 
composite 3-MEKAL model of the diffuse galactic emission 
and the PL a multicolor disk (MCD) component with a very 



Figure 3. Same as Figure l2l for the May 2001 observation (Obs. 
ID 0112290201). 



low temperature (kT — 180 eV; see again Table 0. Adding 
the MCD component to the 3-MEKAL plus PL model leads 
to an improvement in % 2 °f ~ 20. A plot of the [0.8-10.0] 
keV EPIC pn spectrum with the best fitting spectral model 
is shown in Figure [5] 

To check the accuracy of the power-law photon in- 
dex, we tried to fit only the high energy part of the spec- 
trum of M82 X-l in the e nergy interval [3.3-10.0] keV (cfr. 
iFiorito fe Titarchukl20oi) . This should clearly minimize the 
contamination from the galactic diffuse soft emission. In 
this energy range the fits were performed with a power- 
law plus one or two MEKAL components. Without fixing 
the value of the column density, unacceptable fits or unrea- 
sonably high/low values of some parameters are obtained. 
Thus, we kept Nh fixed at the values derived from the best 
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Model 


Column density 
N H 


Temperature 

kT 


Photon index 

r 


Line param. 

energy .width 


Red. X 2 (dof) 






keV 


keV 


Obs.ID 0206080101 


WABS*MEK + 
WABS*MEK + 

WAdo rL 


1.66±0.08 
1.31±0.18 
1.38±0.0 1 


0.51±0.03 
0.99±0.07 


1.54±0.02 




1 Q1 / Q AH\ 

l.ol(o4bj 


WABS*MEK + 
WABS*MEK + 
WABS*(PL+GAUSS) 


1.67±0.08 
1.31±0.19 
1.38±0.07 


0.51±0.04 
0.96±0.07 


1.56±0.02 


E=6.64±0.03 
a =0.09±0.01 


1.27(843) 


WABS*MEK + 
WABS*MEK + 
WABS*MEK + 
WABS*PL 


1.67±0.10 
1.44±0.20 
8.86±4.34 
1.31±0.11 


0.49±0.06 
0.94±0.08 
7.59±1.22 


1.61±0.05 




1.24(843) 


WABS*MEK + 
WABS*MEK + 
WABS*MEK + 
WABS*(PL+MEK) 


1.97±0.34 
1.47±0.16 
9.96±4.64 
1.39±0.10 


0.43±0.06 
0.90±0.03 
8.24±1.21 
0.14±0.03 


1.65±0.05 




1.24(841) 


WABS*MEK + 
WABS*MEK + 
WABS*MEK + 
WABS*(PL+MCD) 


1.60±0.12 
1.68±0.18 
13.72±10.04 
1.76±0.24 


0.48±0.08 
0.92±0.07 
6.35±1.44 
0.18±0.01 


1.66±0.06 




1.22(841) 


Obs.ID 0112290201 




WABS*MEK + 
WABS*MEK + 
WABS*(PL+MCD) 


1.50±0.04 
1.82±0.14 
5.99±0.31 


0.56±0.01 
1.67±0.11 
0.19±0.01 


1.53±0.04 




1.17(884) 



Table 1. Parameters from the spectral fits of M82 X-l in the [0.8 -10-1 keV band. Results are pres ented for our new XMM-Newton 
observation (Obs. ID 0206080101) and for the archival observation of lStrohmaver fc Mushotzkvl i2003fi (Obs. ID 0112290201). The data 
were extracted from a region of 13" radius for both observations (see text for details). 



fits in Table The resulting photon index is in the range 
r = 1.73 — 1.80, in agreement within 2a with that derived 
from the full spectral fit, V = 1.66 (see Tablc0. 

2.2 XMM-Newton observation of May 2001 

The XMM observation taken o n May 6, 2001 (Observa- 
tion ID 0112290201), analyzed bv lStrohmaver fc Mushotzkvl 
(2003), lasted 31 ks. We re-analyzed it following the proce- 
dure outlined in the previous subsection. Standard flares 
filtering leaves 22 and 29 ks of good time intervals for the 
pn and MOS cameras, respectively. Source counts were ex- 
tracted from a 13" region centered on the same source posi- 
tion. The net count rate is ~1.6 count s _1 for pn and ~0.56 
count s _1 for each MOS. The best fit of this observation 
is obtained with a dual MEKAL model plus a PL+MCD 
with different absorbing columns. Adding a third MEKAL 
component leads to an improvement of the fit, but the tem- 
perature converges to unphysically large values (~ 100 keV). 
The parameters of the best fit are reported in Tabled The 
temperature of the MCD component and the PL photon in- 
dex are consistent, within 2a, with those of the April 2004 
observation. However, the absorbing column of these com- 
ponents is approximately 3 times larger than that measured 
in the 2004 observation. It should be noted that the rela- 
tive normalization and temperatures of the MEKAL compo- 



nents vary slightly between the 2001 and 2004 observations. 
Clearly, this is not a consequence of a variation in the dif- 
fuse galactic emission, but reflects simply the fact that the 
model adopted for it is approximate. 

Also in this case, we checked the accuracy of the PL pho- 
ton index fitting the spectra in the restricted energy interval 
[3.3-10.0] keV. The fits were performed with a PL, with or 
without a MEKAL component. As for the 2004 observa- 
tion, without freezing Nh unacceptable fits or unreasonably 
high/low values of some parameters are obtained. Fixing the 
column density at the value reported in Table the result- 
ing photon index is in the range F = 1.66 — 1.68, consistent 
within 3a with that of the full spectral fit, F = 1.53. 

Finally, we note that, although 

IStrohmaver fc Mushotzkvl i2003T) reported evidence for 
a Fe line at 6.55 keV from an analysis of the spectral 
residuals, we are unable to confirm its detection (the 
improvement obtained adding a gaussian component to the 
fit is negligible). 



3 X-RAY LIGHT CURVE 

Table|2]shows the long time scale flux variability of M82 X-l 
from XMM and Chandra data. Measurements from RXTE 
are not reported because it does not have imaging capa- 
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Observatory 


Date 


Count Rate 




Absorbed Flux a 


Unabsorbed flux" 






(count s — 1 ) 




(erg cm -2 s — 1 ) 


(erg cm -2 s —1 ) 


Chandra HRC-I 


1999-10-28 


0.07 6 


7.8 


X 10~ 12 /5.1 X 10~ 12 c 


3.9 x 10 _u /1.4 x 10" 11 c 


Chandra HRC-I 


2000-01-20 


0.52 b 


5. 8 


x 10" n /3.8 x lO" 11 c 


2.9 x 10 -10 /1.1 x 10" 10 c 


XMM EPIC 


2001-05-06 


1.6/0.56 d 




2.4 x 10" 11 


1.2 x 10" 10 


XMM EPIC 


2004-04-21 


1.0/0.35 d 




1.2 x 10" 11 


3.4 x 10" 11 



a For the XMM EPIC observations, average of the pn and MOS fluxes as calculated by XSPEC 
b From Matsumoto et al. C20011 

c See text for details about the adopted spectral model 
d pn/MOS count rates 



Table 2. The 0.2-10 keV fluxes of M82 X-l from the Chandra and XMM observations. 



bilities and contamination from the host galaxy and nearby 
point sources may significantly affect the count rate. For the 
XMM observations the fluxes are those of the best fit mod- 
els of Table Q avaraged over the three EPIC instruments. 
For the Chandra HRC-I observations, the fluxes are calcu- 
lated using the web interface to PIMMS (v. 3.6c). We adopt 
a simplified MEKAL and MCD+PL model with two differ- 
ent absorbing columns and approximate the MEKAL/MCD 
component with a Bremsstrahlung/Blackbody at the same 
temperature. The Chandra fluxes are probably underesti- 
mated because the shape of the MCD and blackbody spectra 
start to differ significantly below 0.5 keV. 

The highest recorded flux emitted by M82 X-l (F ~ 
10 -10 erg cm -2 s _1 ) corres ponds to a luminosity L ~ 2 x 
10 41 erg s" 1 (D = 3.9 Mpc: ISakai fc Madordll999l) . Making 
the usual assumption that, at maximum, the source emits 
at the Eddington limit, we can derive a rough estimate of 
the BH mass, M B h ~ 1500M Q . 



4 TIMING ANALYSIS 

4.1 XMM-Newton data 

For the timing analysis of the 2004 data we avoided the inter- 
val with high background radiation and limited the extrac- 
tion to the longest (nearly) uninterrupted segment of data 
(66 ks) free from solar flares with count rate higher than 
30 count s _1 . To minimize galactic contamination, source 
counts were extracted from a circular region of 8" radius 
and at energies > 2 keV. We produced a light curve from 
pn+MOS data with a time binning of 0.5 s. A few gaps of 
typical duration of ~100 s were present in the light curve and 
were filled with a Poissonian realization around the mean 
value of counts before and after t he gap. We produ ced a 
power spectrum (normalized after iLeahv et ail Il983*l from 
the resulting light curve and rebinned it by a factor of 256 
reaching a frequency resolution of 3.9 mHz (see Figure 
A rather strong QPO peak is evident in the figure. We fit- 
ted the power spectrum with a model consisting of a con- 
stant (for the Poissonian level) plus two Lorentzian com- 
ponents (see lBelloni et ai1l2002t) : one zero-centered for the 
broad band-limited noise and one for the QPO peak. The 
characte ristic frequency for the band-limited noise compo- 
nent (see lBelloni et all2002l for a definition) is 39.4±8.6 mHz 
and its integrated fractional rms is ~22% (after subtracting 
the contribution of the host galaxy). The parameters of the 
QPO can be seen in Table|3| The quality value Q, defined as 



2.8 




Frequency (Hz) 



Figure 4. Total power spectrum from the 2004 XMM observation 
(limited to the range 0.001-1 Hz). The line represents the best fit 
model (see text). 



the ratio of the centroid frequency over the FWHM of the 
QPO, is 4.3±0.5. We repeated the analysis in two separate 
energy bands, 2-4 keV and 4-10 keV. The fractional rms of 
the QPO in these bands resulted to be 13.8% and 23.9% 
respectively. 

In order to investigate the possible variability of the 
QPO during the observation, we produced a spectrogram, by 
aligning power spectra obtained from consecutive stretches 
of data 2048 seconds long. The spectrogram is shown in 
Figure^] For display purposes, in the spectrogram each pixel 
is smoothed as the bilinear interpolation of the values at its 
vertices. A trend towards lower QPO frequencies is apparent, 
correlated with the source count rate (top panel). In order to 
quantify the decrease in centroid frequency, we divided the 
66 ks interval in two segments of 33 ks each and repeated the 
power spectral analysis described above. The two resulting 
power spectra are shown in Figure A fit with the same 
model used for the total power spectrum confirms that the 
centroid frequency of the QPO decreased by 10.8±4.0% (see 
Table 0. 

In Galactic BHCs the frequency of some QPOs is cor- 
related with certain spectral parameters, in particular with 
the power law spectral inde x (see e.g. the QPO fr equency-r 
relation in GRS 1915+105: iVignarca etlITl2003l) . We then 
analyzed the spectra of M82 X-l in the two 33 ks intervals 
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Parameter 


Total observation 


First half 


Second half 


uo (mHz) 


113±2 


120±3 


107±4 


FWHM{mHz) 


26±3 


21±4 


19±3 


Frac % rms 


18.3±1.0 


17.5±1.1 


17.3±1.1 


Signif. (it) 


8.9 


8.3 


8.2 



Table 3. Parameters of the XMM QPO. Errors are at la level. 



Figure 5. Spectrogram obtained from 2048 s XMM data 
stretches. Darker gray corresponds to higher power. [This Fig- 
ure is available as a separate jpg file] . 
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Figure 7. Time history of the centroid frequencies detected from 
M82 X-l in the XMM and RXTE data. The triangles indicate 
the times of RXTE observations when no significant QPO was 
detected. The pairs of dotted lines indicate the range of frequen- 
cies detected by XMM in 2004 and the corresponding intervals at 
half and 1.5 times the frequency. 



Figure 6. (a) Power spectrum from the first half of the XMM 
data, (b) Power spectrum from the second half. The lines show 
the best fit models described in the text. 



after filtering for solar flares, adopting the best fitting model 
reported in Table0 The spectral analysis does not show any 
evidence of variability in the PL spectral index in the two 
intervals. 

We also re-analyzed the XMM o bservation 0112290201 
of IStrohmaver fc Mushotzkvl (120031) . The analysis of the 
power de nsity spectrum confirms the pres ence of the QPO 
found bv IStrohmaver fc Mushotzkvl ll2003l) . 



4.2 RXTE data 

In order to investigate the variability of the QPO frequency 
on longer time scales, we extracted from the RXTE pub- 
lic archive all 30 public observations of M82, spanning over 
the year 1997. For each observation, we accumulated PCA 
light curves in the channel range 0-35, corresponding to 2-13 
keV, with a 0.5 s bin size and produced power spectra in the 
same way as for the XMM data. We detected a significant 
QPO in seven observations , inclu ding t he three reported by 
IStrohmaver fc Mushotzkvl fcOOSt and iFiorito fc Titarchukl 
j2004F ~ These detections are summarized in Table |1| and 



a timing history of their centroid frequencies is shown in 
Figure |7| where also the XMM detections are indicated. Al- 
though the frequencies are variable, they are roughly consis- 
tent with th ree groups in harmonic 1:2: 3 ratio, as recently 
suggested bv lFiorito fc Titarchukl i2004h . In order to calcu- 
late the significance of such an harmonic relation we did a 
numerical simulation and found the nine QPO frequencies 
(the seven from RXTE data plus the two from XMM) to be 
consistent at 2.8 a with being harmonics of a fundamental 
frequency of 54.9 Hz. However, this is not sufficient to com- 
pletely rule out that such a distributions occurs by chance 
and then we did not explore further the consequences of this 
result. More detections are clearly needed in order to address 
this issue. 



5 DISCUSSION 
5.1 Energy spectra 

M82 X-l is one of the best studied ULXs. It suffers from huge 
contamination in the soft X-rays from the diffu se emission of 
the h ost galaxy, the starburst prototype M82 llOriglia et alJ 
120041) . Despite the care in subtracting the contribution of 
this emission, a residual contamination is present in the 
spectra. We modelled this emission with a 3-MEKAL model 
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ObsID 


Date 


Exp. (s) 


QPO v (mHz) 


n a 


20303-02-01-00 


1997 Feb 02 


3709 


166±6 


6.8 


20303-02-02-00 


1997 Feb 24 


3616 


54±5 


5.3 


20303-02-03-00 


1997 May 16 


3312 


50±5 


6.0 


20303-04-04-00 


1997 Jun 07 


3872 


114±5 


5.8 


20303-04-05-00 


1997 Jun 10 


2848 


87±22 


4.4 


20303-08-07-00 


1997 Jul 16 


3127 


67±5 


6.5 


20303-02-04-00 


1997 Jul 21 


2896 


110±2 


4.3 



Table 4. Parameters of the RXTE QPO. Errors are at la level. 



with different absorbing columns. The third MEKAL com- 
ponent at high temperature can account for the He-like Fe 
excess emission in the spectrum, suggesting that it may orig- 
inate from the diffuse emission in M82 and not from the 
ULX. The upper limit to the flux of a line at ~ 6.6 keV in 
the 2001 XMM observation is consistent with the line flux 
measured in the 2004 XMM observation. Since in 2001 the 
source count rate was 50% higher, it might have been dim- 
cult to detect a line with a flux comparable to that of the 
2004 observation. This result appears to be consistent with 
the absence of line flux variability and then with a galac- 
tic origin. However, it should be noted that the very high 
column density of the third MEKAL component may repre- 
sent a problem for this interpretation, although absorption 
columns as high as 10 23 cm" 2 are reported in not-dissimilar 
galaxies iPietsch et al.ll200ll) . 

The HRI instrum ent on board Chandra allowed 
iMatsumoto et all (l200lT) to resolve the central part of the 
galaxy, revealing the presence of a number of point sources 
in addition to M82 X-l (the Chandra source CXOM82 
.10955 50.2+694047, labelled source 7 bv IMatsumoto et alJ 
l200lf) . A total of 9 sources are present in a field of 1' x 1' 
centered on the galaxy. Most of these sources are rather close 
to M82 X-l, below the resolving power of the XMM EPIC 
cameras. For this reason, as mentioned in J2I we subtracted 
from the spectral extraction region circular patches (5") in 
coincidence with the position of the Chandra sources. Al- 
though we are aware that for XMM this is not as effective as 
for Chandra (because a region of 5" corresponds to a XMM 
encircled energy fraction of about 30-50%), this appears to 
be a reasonable compromise between minimizing the con- 
tamination from these sources and preserving the counting 
statistics of M82 X-l. This procedure could not be adopted 
to eliminate the contribu tion from the three near est sources 
to M82 X-l (4, 5 and 6 in lMatsumoto et all2o"oih . In partic- 
ular, source 5 is highly variable. In the Chandra observation 
of October 1999 the count rate of this source was only 30% 
of that of M82 X-l. Later, the former decreased by an or- 
der of magnitude while the latter increased by a factor ~ 7. 
Therefore, while sources 4 and 6 are roughly constant and 
weak with respect to M82 X-l, significant contamination 
from source 5 can not be ruled out. 

A number of ULXs show very soft spectral components 
well fitted by a MCD model, representing emission from an 
accretion disk. Our best fit includes a soft MCD compo- 
nent with temperature kTMCD~0-i8 keV. A plot of the A\ 2 
contours for kT vs. normalisation of the MCD component 
shows that the temperature is well constrained while, within 
3<r, the normalization varies between ~ 500 and ~ 5000 



(XSPEC units; best fit value 1745) . Although this is a large 
variation, the lower bound makes the detection of the MCD 
component quite robust. Adding this component results in 
a statistical improvement of the fit. However, although a 
3-MEKAL model with different absorbing columns offers a 
fair representation of the starbust diffuse emission, it is not 
possible to rule out completely a residual contamination at 
very soft energies from the host galaxy. Under the assump- 
tion that the softest emission is entirely due to an accre- 
tion disk, the temperature of the MCD fit (Tmcd) can be 
used to estimate the BH mass. We note that such an in- 
terpretation in terms of a standard accretion disk does not 
encounter the severe inconsistencies that are found when a 
similar hypothesis is app lied to the soft excess in AGNs (see 
iGierlihski fe Dona |2004). Assuming that Tmcd represents 
an estimate of the maximum temperature of a standard ac- 
cretion dis k, it is (3GM B hM /8 ivar?„) 1/4 = qTmcd, with 
a ~ 2 (see Z ampieri et alJl2004f) . If the disk terminates at 
the innermost stable circular orbit of a Schwarzschild BH, 
the expression for the BH mass then becomes: Mbh/Mq = 
{Mc 2 /L Edd )f 4 {aT M CD/1.5 x 1 7 K)- 4 , where f ~ 1.6 - 
1.7 is a color correction factor jShjmur a fe Takaharalll995l : 
IZampieri. Turolla fe Szuszkiewiczll200 lT) . Assuming Edding- 
ton limited accretion, the BH mass is Mbh ~ 2OO/ 4 M0, 
roughly in agreement with the value inferred from the X- 
ray flux. 



5.2 Timing 

M82 X-l is at present the only ULX where a QPO has been 
discovered. There are a total of 9 detections, two of which 
in the XMM data and seven in the RXTE data, with fre- 
quencies in the range ^50-170 mHz. An important issue is 
of course the possible identification of this QPO with one of 
the QPO types observed in the X-ray light curves of stellar- 
mass BHCs. In the following we will summarize the main 
properties of the QPO in M82 X-l, and discuss its similari- 
ties and differences with the QPOs observed in BHCs. 

• The lowest and highest observed frequencies are 50 ± 5 
mHz and 166 ± 6 respectively (see Figure 0. 

• The frequency distribution over this range is suggestive 
of a harmonic 1:2:3 ratio between them. 

• In the 2004 XMM observation the frequency is observed 
to vary by a factor of ~10% in less than one day. 

• The QPO peak has a quality value higher than 4 (up 
to ~6 in one case). 

• It shows a high fractional rms (up to ~18%). 

• The underlying band limited noise is strong (fractional 
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rms ~22%) and has a characteristic frequency comparable 
to the QPO frequency. 

• The integrated fractional rms of the QPO above 4 keV 
is higher than below that energy. 

Let us compare now all these properties with those of 
the various types of QPOs observed in BHCs. 

• Very-low frequency QPOs. An association of the QPO 
in M82 X-l with the very low frequency (y < 0.02 Hz) 
"heart-beat" QPOs observed in GRS 1915+105 is unlikely 
as their frequency is lower. This, assuming an inverse scaling 
with the black-hole mass, would imply a very low (~solar) 
mass black hole in M82 X-l, which is not in agreement with 
the spectral evidences. 

• High-frequency QPOs. The observed short-time scale 
variability seems to exclude an association with the high- 
frequency "hecto-Hertz" QPOs observed in BHCs, since the 
latter have been detected at rather stable frequencies. How- 
ever, it is worth noticing that, assuming a 10 3 scaling factor 
between the two phenomena, the ~ 10 mHz variation of the 
centroid frequency on a time scale of ~ 30 ks observed in 
M82 X-l would correspond to a ~ 10 Hz variation on a time 
scale of ~ 30 s in a 10 solar mass BHC. Such a short time 
scale variability is impossible to detect in BHCs with the 
present instrumentation. On the other hand, the presence of 
a strong underlying band limited noise, with a characteris- 
tic frequency comparable to the QPO frequency, is clearly 
at variance with the high-frequency "hecto-Hertz" QPOs ob- 
served in BHCs. Furthermore, the rms amplitude of the 
QPO itself in M82 X-l is roughly an order of magnitude 
bigger than that of the "hecto-Hertz" QPOs in BHCs, thus 
making the association very unlikely. For the sake of com- 
pleteness we stress that the detection of the QPO at ~ 166 
mHz reported in this paper lowers the upper limit for the 
mass of the black hole in M82 X-l (assuming that this fre- 
quency is associated with the Keplerian frequency at the 
innermost circular orbit around a Schwarzschild black hole) 
to ~ 1.2 x 10 4 M o . 

• Low-frequency QPOs. Three main types of 
low-frequency QPOs are o bserved in BHCs (see 
ICasella. Belloni fc Stella! 120051) . In two of them, the 
type-A and the type-B QPOs, the peak appears always at 
frequencies near 8 and 6 Hz respectively. Moreover, they are 
both characterized by a weak (a few %) underlying red noise 
component. These properties make an association with the 
variable, strong QPO observed in M82 X-l unlikely. In 
the case of type-A QPOs, its low coherence and amplitude 
make the association even less likely. 

The properties of the QPO in M82 X-l are on the contrary 
reminiscent of those of the third type of BHCs low-frequency 
QPO, the type C, whose characteristic frequencies vary in 
the range 0.1-15 Hz. The similarities in fractional rms, vari- 
ability, quality value, and underlying noise strongly suggest 
an association between the two features. Furthermore, in the 
2004 XMM observation there is evidence for a positive corre- 
lation of the QPO frequency with the count-rate, and a sim- 
ilar correlation is often observed in type-C QPOs. However, 
during the 2001 XMM observation (when the QPO was de- 
tected at a lower frequency) the count rate was higher than 
during the 2004 observation. Since the count-rate vs. fre- 
quency correlation in BHCs is "outburst dependent" (which 
means that during different outburst a source can show simi- 



lar frequencies at different count-rates) no conclusion can be 
derived from the observed phenomenology in M82 X-l. No 
information on the count-rate variability could be obtained 
from the RXTE observations, given the lack of imaging ca- 
pabilities of the satellite. 

Assuming that the QPO detected in M82 X-l is a type- 
C QPO, and scaling the frequency inversely to the BH mass, 
the observed frequency range (from 50 to 166 mHz) would 
yield Mbh anywhere in the range 10-1000 Mq. However, 
type-C QPOs are observed in BHCs throughout the whole 
Hard-Intermediate State fsee lHoman fc Bellonil 12005). and 
their frequency is known to decrease with the hardness of 
the energy spectrum. At the lowest observed frequencies, 
the spectrum is hard and there is often no evidence for the 
presence of a soft thermal component. As the contribution 
from a disk appears and increases, the QPO frequency also 
increases. The XMM spectra of both observations in which a 
QPO has been detected in M82 X-l show possible evidence 
for a disk contribution (see before). To the extent that the 
two phenomena can be compared, the presence of a soft 
component would exclude that the type-C QPO in M82 X-l 
is in the lowest frequency range, increasing the lower limit 
for Mbh. 



6 CONCLUSIONS 

In this paper we reported a complete analysis of XMM- 
Newton and RXTE observations of the superluminal X-ray 
source M82 X-l. The similarities in fractional rms, vari- 
ability, quality value, and underlying noise strongly suggest 
an association between the QPO in M82 X-l and the low- 
frequency, type-C QPOs observed in BHCs. The combined 
spectral and timing analysis allows us for the first time to 
put strong constrains to the mass of the central black hole 
in this source, yielding to a value between a few tens to one 
thousand solar masses. 
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